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Abstract 

In this paper, we propose an evolving Sierpinski gasket, based on which we estabUsh 
a model of evolutionary Sierpinski networks (ESNs) that unifies deterministic Sier- 
pinski network [Eur. Phys. J. B 60, 259 (2007)] and random Sierpinski network [Eur. 
Phys. J. B 65, 141 (2008)] to the same framework. We suggest an iterative algo- 
rithm generating the ESNs. On the basis of the algorithm, some relevant properties 
of presented networks are calculated or predicted analytically. Analytical solution 
shows that the networks under consideration follow a power-law degree distribution, 
with the distribution exponent continuously tuned in a wide range. The obtained 
accurate expression of clustering coefficient, together with the prediction of average 
path length reveals that the ESNs possess small- world effect. All our theoretical 
results are successfully contrasted by numerical simulations. Moreover, the evolu- 
tionary prisoner's dilemma game is also studied on some limitations of the ESNs, 
i.e., deterministic Sierpinski network and random Sierpinski network. 
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Fig. 1. The first two stages of construction of the Sierpinski gasket (a) and its 
corresponding network (b). 



1 Introduction 



In the last few years, complex networks have attracted a growing interest 
from a wide circle of researchers [T|2|3P] . The reason for this boom is that 
complex networks describe various systems in nature and society, such as the 
World Wide Web (WWW), the Internet, collaboration networks, and sexual 
network, and so on. Extensive empirical studies have revealed that real-life 
systems have in common at least two striking statistical properties: power- 
law degree distribution [5], small- world effect [6] including small average path 
length (APL) and high clustering coefficient. In order to mimic real- word sys- 
tems with above mentioned common characteristics, a wide variety of models 
have been proposed p]l2|3|^ . At present, it is still an active direction to con- 
struct models reproducing the structure and statistical characteristics of real 
systems. 

In our previous papers, on the basis of the well-known Sierpinski fractal (or 
Sierpinski gasket), we have proposed a deterministic network called determin- 
istic Sierpinski network (DSN) [7J, and a stochastic network named random 
Sierpinski network (RSN) [8J, respectively. Both the DSN and RSN possess 
good topological properties observed in some real systems. In this paper, we 
suggest a general scenario for constructing evolutionary Sierpinski networks 
(ESNs) controlled by a parameter q. The ESNs can also result from Sierpinski 
gasket and unify the DSN and RSN to the same framework, i.e., the DSN 
and RSN are special cases of RSNs. The ESNs have a power-law degree distri- 
bution, a very large clustering coefficient, and a small intervertex separation. 
The degree exponent of ESNs is changeable between 2 and 3. Moreover, we 
introduce a generating algorithm for the ESNs which can realize the construc- 
tion of our networks. In the end, the cooperation behavior of the evolutionary 
prisoner's dilemma game on two limiting cases (i.e., DSN and RSN) of the 
ESNs is discussed. 
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Fig. 2. (Color online) The sketch maps for the construction of random Sierpinski 
gasket (left panel) and its corresponding network (right panel). 



2 Brief introduction to deterministic and random Sierpinski net- 
works 



We first introduce Sierpinski gasket, which is also known as Sierpinski triangle. 
The classical Sierpinski gasket denoted as St after t generations, is constructed 
as follows piTO] : start with an equilateral triangle, and denote this initial 
configuration as 5*0. Perform a bisection of the sides forming four small copies 
of the original triangle, and remove the interior triangles to get Si. Repeat 
this procedure recursively in the three remaining copies to obtain 5*2, see 
Fig. [D^a). In the infinite t limit, we obtain the famous Sierpinski gasket St- 
From Sierpinski gasket we can easily construct a network, called deterministic 
Sierpinski network, with sides of the removed triangles mapped to nodes and 
contact to edges between nodes [7]. For uniformity, the three sides of the 
initial equilateral triangle at step also correspond to three different nodes. 
Figure [T](b) shows a network based on 5*2. 

Analogously, one can construct the random Sierpinski network [8] derived 
from the stochastic Sierpinski gasket, which is a random variant of the deter- 
ministic Sierpinski gasket. The initial configuration of the random Sierpinski 
gasket is the same as the deterministic Sierpinski triangle. Then in each of the 
subsequent generations, an equilateral triangle is chosen randomly, for which 
bisection and removal are performed to form three small copies of it. The 
sketch map for the random fractal is shown in the left panel of Fig. [21 From 
this fractal we can easily establish the random Sierpinski network with sides 
of the removed triangles mapped to nodes and contact to links between nodes. 
The right panel of Fig. [2] gives a network derived from the random Sierpinski 
gasket. 
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Fig. 3. (Color online) Iterative construction method for the network. 
3 Unifying model and its iterative algorithm 

In this section, we introduce an evolving unified model for the deterministic 
and random Sierpinski networks. First we give a new variation, called evolving 
Sierpinski gasket (ESG), for the Sierpinski gasket. The initial configuration of 
the ESG is the same as the deterministic Sierpinski gasket. Then in each of 
the subsequent generations, for each equilateral triangle, with probability q, 
bisection and removal are performed to form three small copies of it. In the 
infinite generation limit, the ESG is obtained. In a special case g = 1, the 
ESG is reduced to the classic deterministic Sierpinski gasket. If q approaches 
but is not equal to 0, it coincides with the random Sierpinski gasket described 
in Ref. [8]. The proposed unified model is derived from this ESG: nodes rep- 
resent the sides of the removed triangles and edges correspond to contact 
relationship. As in the construction of the deterministic and random Sierpin- 
ski networks [7l8], the three sides of the initial equilateral triangle (at step 0) 
of the ESG are also mapped to three different nodes. 

In the construction process of the ESG, for each equilateral triangle at ar- 
bitrary generation, once we perform a bisection of its sides and remove the 
central down pointing triangle, three copies of it are formed. When build- 
ing the unifying network model, it is equivalent that for each group of three 
newly-added nodes, three new triangles are generated, which may create new 
nodes in the subsequent generations. According to this, we can introduce an 
iterative algorithm to create the ESNs. Using the proposed algorithm one can 
write a computer program conveniently to simulate the networks and study 
their properties. 

We denote the ESNs after t (t > 0) iterations by EW(t), then the proposed 
algorithm to create ESNs is as follows. Initially {t = 0), EW{0) has three 
nodes forming a triangle. At step 1, with probability q, we add three nodes 
into the original triangle. These three new nodes are connected to one another 
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shaping a new triangle, and both ends of each edge of the new triangle are 
linked to a node of the original triangle. Thus we obtain EW{1), see Fig. [31 For 
t > 1, EW{t) is obtained from EW{t — 1). For the convenience of description, 
we give the following definition: for each existing triangle in EW{t — l), if there 
is no node in its interior and among its three nodes there is only one youngest 
node (i.e., the other two are strictly elder than it), we call it an active triangle 
(with the initial triangle as an exception). At step t — 1, for each existing 
active triangle, with probability q it is replaced by the connected cluster on 
the right of Fig. [31 then EW{t) is produced. The growing process is repeated 
until the network reaches a desired order. When q = 1, the network is exactly 
the same as the DSN [7]. If g < 1, the network grows randomly. Especially, 
as q approaches zero and does not equal zero, the network is reduced to the 
RSN studied in detail in Ref. [8]. 

Next we compute the order (number of all nodes) and size (number of all 
edges) of EW{t). Denote L^it) as the number of active triangles at step t. 
Then, La(0) = 1. By construction, we can easily derive that L^{t) = (H-2g)*. 
Let Ly{t) and Le(t) be the number of nodes and edges created at step t, 
respectively. Note that each active triangle in EW{t — 1) will (see Fig. [3]) lead 
to three new nodes and nine new edges in EW{t) with probability q. Then, at 
step 1, we add expected L^i^l) = 3q new nodes and Le{l) = 9q new edges to 
EW{0). After simple calculations, one can obtain that at step {ti > 1) the 
number of newly-born nodes and edges is L^iti) = 3q L^(ti — 1) = 3q{l+2qy^~^ 
and Le(ti) = 9qL^(ti — 1) = 9g(l + 2g)*'~^, respectively. Thus the average 
number of total nodes Nt and edges Et present at step t is 

ti=l ^ 

and 

^..3.i:M,)^?(i±Mi^. (2) 

ti=\ ^ 

respectively. So for large t, the average degree kt = ^ is approximately 6. 
Obviously, we have Et = 3Nt — 6. Moreover, according to the connection rule, 
arbitrary two edges in the ESNs never cross each other. Thus, the class of 
networks under consideration are maximal planar networks (or graphs) 



4 Relevant characteristics of the Networks 

In the following we will study the topology properties of EW{t), in terms of 
degree distribution, clustering coefficient, and average path length. 
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4-1 Degree distribution 



When a new node i is added to the network at step tj {ti > 1), it has a 
degree of 4. Let Li^{i,t) be the expected number of active triangles at step 
t that will create new nodes connected to the node i at step t + 1. Then at 
step ti, L^{i,ti) = 1. From the iterative generation process of the network, 
one can see that at any subsequent step each two new neighbors of i generate 
two new active triangles involving i, and one of its existing active triangles is 
deactivated simultaneously. We define ki(t) as the degree of node i at time t, 
then the relation between ki{t) and Li^{i,t) satisfies: 

^.(M) = ^^^. (3) 



Now we compute Li^{i,t). By construction, L^^ii^t) = (1 + q)L^{i,t — 1). 
Considering the initial condition L^{i,ti) = 1, we can derive Li^{i,t) = (1 + 
qy-U _ Then at time t, the degree of vertex i becomes 

A;,(t) = 2(l + g)*-*'+2. (4) 

Since the degree of each node has been obtained explicitly as in Eq. (j4]), 
we can get the degree distribution via its cumulative distribution [3], i.e., 
-Pcum(^) = Y.k'>k^i.k' I't) l^t ~ k^~'^ , where N{k',t) denotes the number of 
nodes with degree k'. The detailed analysis is given as follows. For a degree 
k = 2(1 + g)*-™ + 2, there are Ly{m) = 3g(l + 2q)'^~^ nodes with this exact 
degree, all of which were born at step m. All nodes born at time m or earlier 
have this or a higher degree. So we have 

E Nik', t)=A + tus) = + + (5) 

k'>k s=l 

Thus, the cumulative degree distribution is give by 

PcUk) = ^ E N{k', t) = (6) 
Nt 3(1 + 2g)* + 3 

In^ 

Substituting for m in this expression using m = t — i^^^il^^ gives 

3(l + 2g)*-2W+ir(fc-2) W+ir+3 
3(l + 2g)* + 3 ■ 

When t is large enough, one can obtain 

lii(l+2.;) ln(l+29) 

Pcum(^) ~2M^(A;-2) (8) 
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So the degree distribution follows a power-law form P{k) k with the 
exponent 7 = 1 + ^j^^^^- Note that the degree exponent 7 is a continuous 

function of q, and belongs to the interval (1 + 3]. As g decreases from 1 to 
0, 7 increases from 1 + |^ to 3. In the two limitations, i.e., q = 1 and g — 
(but g 7^ 0), the evolutionary Sierpinski networks reduce to the deterministic 
Sierpinski network [7j and its stochastic variant [H], respectively. Figure H] 
shows, on a logarithmic scale, the scaling behavior of the cumulative degree 
distribution Pcum(^) for different values of q. Numerical simulation agrees well 
with the analytical result. 




Fig. 4. (Color online) The cumulative degree distribution Pcum{k) at various q values. 
The evolutionary steps of simulated network for q = 0.8, q = 0.3, and q = 0.001, 
are t = 12, t = 26, and t = 200000, respectively. The three lines are the theoretical 
results as provided by Eq. ([8]). All data are from the average of ten independent 
simulations. 



4-2 Clustering coefficient 



In a network, the clustering coefficient j6j Ci of node i is defined as the ratio 
between the number of edges Cj that actually exist among the fcj neighbors of 
node i and its maximum possible value ki{ki — l)/2, i.e., Ct = 2ei/[ki{ki — 1)]. 
The clustering coefficient of the whole network is the average of C-s over all 
nodes in the network. 

For the ESNs, the analytical expression of clustering coefficient C{k) for a 
single node with degree k can be derived exactly. When a node is added into 
the network, its ki and Cj are both 4. At each subsequent discrete time step, 
each of its active triangles increases both ki and Cj by 2 and 3, respectively. 
Thus, Cj = 4 + |(A;j — 4) for all nodes at all steps. So there is a one-to-one 
correspondence between the clustering coefficient of a node and its degree. For 
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Fig. 5. The clustering coefficient of the whole network as a function of the size of 
the network for various q. Results are averaged over ten network realizations for 
each datum. 

a node of degree /c, we have 



C{k) 



4 + f(A;-4) 



k{k-l) 



4 



k k-r 



(9) 



which is inversely proportional to k in the limit of large k. The scaling of 
C{k) ~ k~^ has been empirically observed in many real-life networks [12] . 

After t generation evolutions, the average clustering coefficient Ct of network 
EW{t) is given by 



1 ' 

— E 

Nth 



4 



1 



LJr) 



(10) 



where the sum runs over all the nodes and is the degree of those nodes 
created at step r, which is given by Eq. (jlj). In the infinite network order limit 
{Nt — > oo), Eq. (ITU]) converges to a nonzero value C, see Fig. O Moreover, it 
can be easily proved that both Ct and C increase with q. Exactly analytical 
computation shows: when q increases from to 1, C grows from 0.5674 [8] to 
0.5745 [7J. Therefore, the evolutionary networks are highly clustered. Figure 
E] shows the average clustering coefficient of the network as a function of q, 
which is in accordance with our above conclusions. From Figs. S] and [HI one 
can see that both degree exponent 7 and clustering coefficient C depend on 
the parameter q. The mechanism resulting in this relation deserves further 
study. The fact that a biased choice of the active triangles at each iteration 
may be a possible explanation, see Ref. 
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Fig. 6. The dependence of the average clustering coefficient on parameter 
4-3 Average path length 



From above discussions, one knows that the existing model shows both the 
scale-free nature and the high clustering at the same time. In fact, our model 
also possesses small- world property. Next, we will show that our networks have 
at most a logarithmic average path length (APL) with the number of nodes, 
where APL means the minimum number of edges connecting a pair of nodes, 
averaged over all couples of nodes. 

Using a mean-field approach similar to that presented in Refs. P^ffTB] , one can 
predict the APL of our networks analytically. By construction, at each time 
step, the number of newly-created nodes is different. In order to distinguish 
different nodes, we construct a node sequence in the following way: when AA^ 
nodes are added at a given time step, we label them as M + 1, M + 2, . . . , M + 
AN, where M is the total number of the pre-existing nodes. Eventually, every 
node is labeled by a unique integer, and the total number of nodes is Nt = 
3(i+2^g)'+3 ^j^g ^_ denote L{N) as the APL of the ESNs with order 

A^. It follows that L{N) = where D{N) = J2i<i<j<N dij is the total 

distance, and where dij is the smallest distance between node i and node j. 
Note that the distances between existing node pairs are not affected by the 
addition of new nodes. As in the analysis of p^|IT5] . we can easily derive that 
D{N) ~ N"^ in the infinite limit of A^. Then, L{N) ~ InA^. Thus, there is a 
slow growth of the APL with the network order A^. This logarithmic scaling 
of L{N) with network order A^, together with the large clustering coefficient 
obtained in the preceding subsection, shows that the considered graphs have 
a small- world effect. 



In particular, in the case of g = 1, we can exactly compute the average path 
length. A previously reported analj^ical result has shown that the APL for 
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this special case grows logarithmically with the order of the network [7]. In 
Fig. [TJ we report the simulation results on the APL of ESNs for different q. 
From Fig. [71 one can see that APL decreases with increasing q. For all g, APL 
increases logarithmically with network order. 
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Fig. 7. Semilogarithmic graph of the APL versus the network order N . 

5 Prisoner's Dilemma game on two limiting cases for ESN 

The ultimate goal of study of network structure is to study and understand 
the workings of systems built upon those networks [3|l4] . Recently, some re- 
searchers have focused on the analysis of functional or dynamical aspects of 
processes occurring on networks. One particular issue attracting much atten- 
tion is using evolutionary game theory to analyze the evolution of cooperation 
on different types of networks [16]. Cooperation is ubiquitous in the real-life 
systems, ranging from biological systems to economic and social systems ^7j. 
Understanding the emergence and survival of cooperative behavior in these 
systems has become a fundamental and central issue. 

After studying the relevant characteristics of network structure, which is de- 
scribed in the previous section, we will study the evolutionary game behavior 
on the networks, with focus on the game of prisoner's dilemma (PD). In the 
simple, one-shot PD game, both receive R under mutual cooperation and P 
under mutual defection, while a defector exploiting a cooperator gets amount 
T and the exploited cooperator receives S, such that T > i? > P > 5*. As a re- 
sult, it is better to defect regardless of the opponent's decision, which in turns 
makes cooperators unable to resist invasion by defectors, and the defection is 
the only evolutionary stable strategy in fully mixed populations. 

We now investigate the evolutionary PD game on our networks to reveal the 
influences of topological properties on cooperation behavior. Here we only 
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study two limiting cases: q = I and q (but g 7^ 0). As usual in the 
studies [HUn], we choose the payoffs a.s R = 1, P = S = 0, and T = b > 1, 
and implement the finite population analogue of replicator dynamics. During 
each generation, each individual i plays the single given game with all its 
neighbors, and their accumulated payoff being stored in Pj. After each round 
of the game, the individual i is allowed to update its strategy by selecting at 
random a neighbor among all its neighbors, j, and comparing their respective 
payoffs Pi and Pj. If Pj > Pj , the individual i will keep the same strategy for 
the next generation. On the contrary, the individual i will adopt the strategy 
of its neighbor j with a probability dependent on the payoff difference {Pj — 
Pi) as Hi^j = {Pj — Pi)/ maxjfcj, fcjjfe. All individuals update its strategies 
synchronously during the evolution process. 




1.0 1.2 1.4 1.6 1.8 2.0 

Temptation to defect b 

Fig. 8. Frequency of cooperators as a function of the temptation to defect b. 

In Fig. [HI we report the simulated results (i.e. the dependence of the equilib- 
rium frequency of cooperators on the temptation to defect b) for both DSN and 
RSN. Simulations are performed for both networks with order 9843. Each data 
point is obtained by averaging over 100 simulations for each of ten different 
network realizations. From Fig. [HI one can see that for b < 1.6 the cooperators 
in the deterministic Sierpinski network is dominant over defectors. And simi- 
lar phenomenon is also observed for its corresponding random version. Thus, 
both of the network structures are in favor of cooperation upon defection for 
a wide range of b. Figure [8] also shows that in both networks the frequency of 
cooperators makes a steady decrease when the temptation changes from 1.6 
to 1.75, and then drops dramatically when the temptation increases to 2. On 
the other hand, for large b (such as 6 > 1.6), the equilibrium frequency of DSN 
is higher than that in RSN, and this phenomenon is more obvious when the 
temptation is larger and goes up to 2. 

The observed phenomena in Fig. [8] can be explained according to the underly- 
ing network structures. Since both DSN and RSN are scale-free networks, this 
heterogeneous network architecture makes cooperation become the dominat- 
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ing trait over a wide range of temptation to defect b [T91I20] . Although both 
networks have scale-free property, DSN is more heterogeneous than RSN since 
the former has a smaller exponent of power-law degree distribution than the 
latter; at the same time, the average clustering coefficient of DSN is larger than 
that of RSN. These two different characteristics between DSN and RSN can 
account for the dissimilar cooperation behavior in both networks p^|2T] : A 
higher value of average clustering coefficient, together with a smaller exponent 
of power-law degree distribution, produces an overall improvement of cooper- 
ation in DSN, even for a very large temptation to defect, which is compared 
to that in RSN. 



6 Conclusions 

In summary, on the basis of Sierpinski gasket, we have proposed and studied 
one kind of evolving network: evolutionary Sierpinski networks (ESNs). Ac- 
cording to the network construction process we have presented an algorithm 
to generate the networks, based on which we have obtained the analytical 
and numerical results for degree distribution, clustering coefficient, as well as 
average path length, which agree well with a large amount of real observa- 
tions. The degree exponent can be adjusted continuously between 1 + |^ and 
3, and the clustering coefficient is very large. Moreover, we have studied the 
evolutionary PD game on two limiting cases of the ESN. 

It should be stressed that the network representation introduced here is con- 
venient for studying the complexity of some real systems and may have wider 
applicability. For instance, a similar recipe has been recently adopted for in- 
vestigating the navigational complexity of cities [22]; on the other hand, it is 
frequently used in RNA folding research [23||2^ : moreover, earlier links asso- 
ciating this network representation with polymers have proven useful to the 
study of polymer physics [25)126] . Thus, our study provides a paradigm of rep- 
resentation for the complexity of many real-life systems, making it possible 
to study the complexity of these systems within the framework of network 
theory. 

Because of its three important properties: power-law degree distribution, small 
intervertex separation, and large clustering coefficient, the proposed networks 
possess good structural features in accordance with a variety of real-life net- 
works. Additionally, our networks are maximal planar graphs, which may be 
helpful for designing printed circuits [TT]. Finally, it should be mentioned that 
although our model can reproduce a few topological characteristics of real-life 
systems, it remains unknown whether the model can capture a true under- 
lying mechanism responsible for those properties observed in real networks. 
This belongs to the issue of model evaluation, which is beyond the scope of 
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the present paper but deserves further study in future [27]. 
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